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Abstract: We present a three-body formalism describing the final-state interaction effects 
in the pp — ► ppr\ reaction close to threshold. We derive a three-body enhancement factor 
devised in such a way that all three pair-wise interactions are regarded on equal footing. The 
enhancement factor is obtained by expanding the three-particle wave function in hyperspherical 
harmonics. It has been shown that close to threshold the p — p interaction strongly dominates 
whereas the r\ — p interaction gives almost negligible contribution to the calculated effective 
mass spectra. Within the presented three-body approach it has been possible to reproduce 
the effective mass distributions at the excitation energy Q=15.5 MeV in good accord with the 
£>. ■ data. 
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O ! 1 Introduction 

The last decade has seen major advances in the experimental investigation of the near threshold 
meson production reactions in nucleon-nucleon collisions (for a review cf. 1 and 2 ). In 
the recent measurements of the pp — > ppr) reaction a very accurate determination of the four- 
momenta of both outgoing protons allowed for the full reconstruction of the kinematics of the 
final r\pp state. In consequence, these measurements provided in addition to the r\ and the proton 
angular distributions, also the pp and r/p effective mass distributions [Sill]- The common feature 
of the near-threshold meson production in proton-proton collisions is the dominance of the very 
strong proton-proton final state interaction (FSI). This effect is clearly visible in the invariant 
mass distributions: as a prominent peak close to threshold in the (pp)-mass distribution, or as a 
bump near the end-point in the (rip)-m&ss spectrum. For sufficiently low excitation energies the 
3 Po -^ 1 So s transition amplitude becomes necessarily the sole contributor to the cross section as 
it is the only amplitude surviving at threshold. The supposition that this happens at the lowest 
available excitation energy equal Q=15.5 MeV [I] appears to be quite plausible, especially that 
in a similar experiment [H] at Q=15 MeV the measured angular distributions were consistent 
with isotropy. The description in terms of a simple model in which a constant r\ production 
amplitude is multiplied by pp FSI enhancement factor, although qualitatively correct, is not fully 
satisfactory in quantitative terms. The calculated invariant mass distributions are presented in 
Fig- HI (dashed curves). In order to improve the agreement with experiment two possibilities 
have been considered. The most obvious explanation admits the contribution from the p-waves, 
or, more precisely, from the 5o — ► 3 Pos and 1 Z?2 — ► 3 P2S amplitudes. These amplitudes have the 
best chance to show up when the relative momenta of the final-state protons take the largest 
values allowed by the phase space. Since at Q=15.5 MeV this sector still overlaps with the 
peak region of the ^o enhancement factor, the s-wave also receives there maximal amplification. 
Therefore, the relative strength of the p-wave amplitudes to be discernible has to be quite 
substantial which in general should be reflected by a pronounced angular dependence of the 

1} Talk given at COSY-11 meeting, Cracow, 1-3 June, 2004 
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Figure 1: Invariant mass distributions at Q=15.5 MeV; left panel: (pp)-mass plot; right panel: (rrp)- 
mass plot. The data are from 0] the calculated curves are from |SJ. 



cross section. This difficulty has been thoroughly examined by Nakayama et al. [S] who pointed 
out that the unwanted angular dependence might still be suppressed under two circumstances: 
(i) if the 1 Z?2 — ► 3 P2S amplitude was negligibly small so that the angular dependent term was 
absent, or, (ii) if the lack of angular dependence resulted from cancellations - from a destructive 
interference between ^o — ► 3 Pos and l D2 — » 3 P2S amplitudes. Thus, a model basing upon a 
strong p-wave needs additionally somewhat fortuitous coincidences. According to the second 
proposition presented in |jj|, the s-wave amplitude dominates because it is proportional to the 
huge FSI enhancement factor whilst the p-wave amplitudes are not. Therefore, the latter might 
be neglected explaining in a natural way the lack of angular dependence. Instead, a weak energy 
dependence is admitted in the production amplitude. Both models are capable of improving the 
agreement with experiment at the expense of introducing an adjustable parameter. The best fit 
from jHj is depicted in Fig. ^by the full line. Since the agreement with experiment presented in 
[2] is of similar quality, polarization experiments are required to discriminate between these two 
models. The prediction of the s-wave model [5] is that all polarization observables are bound to 
vanish. Any non-zero value for the analyzing power reveals the presence of higher partial waves 

0. ' 

In the above considerations final-state r/p interaction has been ignored and an interesting 
question arises how much its inclusion can change the resulting invariant mass spectra. The r/p 
interaction is poorly known but there have been suggestions that the corresponding scattering 
length might be as large as about 1 fm. This is still one order of magnitude smaller than the pp 
scattering length but the discrepancies in the invariant mass distribution in Fig. n are n °t large 
either. The purpose of this work is to shed some light on the possible role of the r]p final-state 
interaction but to be able to do that we need a formalism in which all pair-wise interactions 
would be regarded on equal footing. 

The plan of our presentation is as follows. In the next Section we start with the well known 
two-body case recalling the arguments leading to the derivation of the FSI enhancement factor. 
In Section 3 we generalize these ideas to the three-body case by utilizing the hypersherical 



harmonics approach. Finally, in Section 4 we verify our model by confronting the obtained 
results with the experimental data. 

2 Two-body final-state interaction 

Since the proton-proton final-state interaction is believed to be the dominant ingredient in the 
description of the pp — ► ppr] reaction close to threshold, it is logical to begin with the two-body 
FSI problem. The basic idea how to account for final state interaction was put forward 50 years 
ago by Fermi, Watson, Migdal 7 and others (for a review cf. [S]) and is based on the observation 
that in many processes the interaction responsible for carrying the system from the initial to 
the final state is of such a short range that in the first approximation may be regarded as point 
like. As a prototype one may consider a meson (x) production reaction NN — > NNx. To 
generate the meson mass m in nucleon-nucleon collision a large momentum transfer is required 
between the initial and the final nucleons, which is typically of the order \/Mm, with M being 
the nucleon mass. The corresponding "range" of the production interaction is therefore much 
shorter than the range of the interaction between the two final state nucleons. Although it is 
perfectly true that the final state NN interaction significantly distorts the NN wave function but 
in the transition matrix element the contribution from all but the smallest NN separations will 
be strongly suppressed and the main effect may be attributed to the change of the normalization 
of the wave function at zero separation. If the non-interacting NN pair is described by a plane 
wave e~ , where k is the relative NN momentum (h = c = 1 units are used hereafter), to 
account for final state interaction the latter must be replaced in the transition matrix element by 
the complete NN wave function ty~(k, ry satisfying outgoing spherical wave boundary condition 
at infinity. Nevertheless, for a point-like interaction, we may set 

*"(fc,r) t &e- lk - r C(k) (1) 

in the matrix element so that the final state interaction will be accounted for by multiplying the 
transition matrix element by the enhancement factor, defined as 

C(k) = lim * + (-fc,r)/e- jfc - r . (2) 

The factor |C(/c)| 2 that appears in the cross section represents the ratio of two probabilities: one 
of finding the interacting NN pair at zero separation, while the other probability is associated 
with non-interacting particles. By construction, when the final state interaction is turned off, the 
enhancement factor will be equal to unity. Expanding both, the numerator and the denominator 
on the right hand side of (j2j) in partial waves, we have 

C(k) = lim TZ^i + ^'i,Ak,r)lrP, {k .f) 
~° E£o(« + l)'-*J/(*r)^(*-*') 

where ip£(k,r) ~ r i+1 for small r, jc(kr) is spherical Bessel function and Pg(k ■ r) denotes Leg- 
endre polynomial. Clearly, in the limit r — > in (J3J), all higher partial waves will be suppressed 
by the centrifugal barrier, and only the contribution the from s-wave survives. Thus, we obtain 
a simple formula 

C(k) = ifo(k,0)'/k, (4) 

where prime denotes derivative with respect to r, and, as apparent from Ijljl. the enhancement 
factor is determined by the slope of the wave function at the origin. To find this slope we must 
know the NN s-wave interaction and for simplicity we shall in the following assume that the 



latter takes the form of a spherically symmetric radial potential. The shape of this potential 
may be arbitrary but it must be of a short range. Given the NN potential, we can integrate 
outward the appropriate wave equation, containing both the nuclear and the Coulomb potential, 
generating numerically a regular solution uo(k, r) (i.e. vanishing at the origin) whose derivative 
for later convenience is selected to be 

u (k,0y = C ( V )k, (5) 

where n denotes the Sommerfeld parameter and Cq{t]) 2 = 2-irr]/[ex.p(2Trr]) — 1] is the Coulomb 
barrier penetration factor. The sought for physical solution ipo(k,r) occurring in (J3J which is 
also regular, is necessarily proportional to uo(k,r), and, more explicitly, we have 

Mk,r) = [C(k)/C ( V )]u (k,r). (6) 

Now, all we need to calculate C(k) is the asymptotic expression for the physical wave function. 
For r = R with R much bigger than the range of the nuclear potential, the physical wave function 
takes the form 

Mk,R) = F (v,kR) + f (k) H+( v ,kR), (7) 

where Hq (rj, kR) = Go(j],kR) + iFo(n,kR) with Go(rj,kR) and Fo(r),kR) being the standard 
Coulomb wave functions defined in [5j, and fo(k) = sinJe denotes the s-wave scattering am- 
plitude with 5 being the s-wave Coulomb distorted phase shift. The differentiation of ® with 
respect to R, provides us with a second condition for the derivatives but it should be noted that 
uo(k,R) and u(k,R)' occurring in these two matching conditions are to be regarded as known 
quantities. Indeed, they are fully specified by the boundary condition at the origin © and can 
be either calculated analytically, or obtained by numerical methods. Therefore, we end up with 
two algebraic equations in which the two unknowns are the enhancement factor C{k) and the 
scattering amplitude fo(k) and the respective solutions, can be conveniently written as 

C(/C) = w[H+(r ] ,kR),u (k,R)Y (8) 

and 

, (k) = w[F (n,kR),u (k,R)} 

M ' w[H+( V ,kR),u (k,R)V 

where the symbol w[f, g] denotes the Wronskian defined as w[f, g] = fg' — f'g. The specific 
Wronskian occurring in the denominators of (JSJ (EJ has been referred to as the Jost function 
jH]- Thus, given the potential, the two-body enhancement factor can be obtained from (jHJ). 

3 Three-body final-state interaction 

We wish now to extend the ideas outlined above to the three-body case. We assume that the 
pair-wise short range interaction between the particles is strong but the system is Borromean, i.e. 
no binary bound state may exist in the three-body system under consideration. A formal theory 
of the continuum in a Borromean system was developed in jlUj using the same hyper spherical 
harmonics method which has been widely employed for the investigation of bound states and 
specifically the halo nuclei jllj . In this paper we follow this approach to derive the enhancement 
factor for three interacting particles. The wave functions in the continuum are solutions to 
the three-body problem satisfying the correct boundary conditions at infinity where the three- 
body asymptotics is most naturally expressed in terms of the rotationally and permutationally 
invariant hyperradius p defined as the square root of the sum of squares of the inter-particle 



distances. Since the hyperradius reflects the size of the three-body system, similarly as in the 
two-body case, the enhancement factor may be obtained as the limiting value p — > of the three- 
body wave function. The method of hyperspherical harmonics has been well documented in the 
literature, but to make this paper self-contained we wish to summarize briefly the theoretical 
framework necessary to treat the continuum of a Borromean system. 

3.1 Coordinate sets and hyperspherical harmonics 

For assigned particle positions (ri, r2, i^) and masses (mi,m2, "23), the translationally invariant 
normalized sets of Jacobi coordinates Xj,yj are defined, as follows 



Xi = \ 1 7 ^( r j- r k), (10a) 

V (rrij + m k ) 11 



/ m, (rrij + m k ) ( rrij Yj + m k r k \ 

y fiM \ rrij + m k J 

R = (miTi + nijYj + m k r k )/M, (10c) 

where {i,j,k} is a permutation of the particle labels {1,2,3}, M = m\ + rri2 + 7713, and \x 
denotes an arbitrary mass which just sets the mass scale. Each of the three equivalent pairs 
(xj, y$) together with the center of mass coordinate R describes the system. The transformation 
between different sets of Jacobi coordinates is referred to as a kinematical rotation and takes 
the form 

' X A == f-cosu>ij smuij\ /xA 
yj \-smuij - cos tvij J \jj 

where the rotation angle, confined by —tt/2 < tOij < vr/2, is given by 



arctan 



a{i, j, k}^M m k /mi rrij 



(12) 



with cr{i,j, k} denoting the sign of the permutation {i,j, k}. 

The Jacobi momenta kj,qj and P, canonically conjugate to Xj,yj and R, respectively, are 
defined by the relations, 

/ rn ■ rn 1 it l x\ • x\ i \ 

(13a) 
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(13b) 

(13c) 

where p« (i = 1,2,3) are the laboratory frame momenta. 

Instead of the Jacobi coordinates we shall use hyperspherical coordinates comprising of a 
hyperradius p and five angles. The hyperradius p determines the size of a three-body system 
and is invariant with respect to kinematic rotations 

p 2 = x 2 l +yl 1 = 1,2,3. (14) 

The five angular variables forming a five-dimensional solid angle Q p include the usual angles 
(9 l x , (j) l x ), (6y, 4>y) specifying the unit vectors Xj, y^ and these are supplemented by the hyperangle 
at defined by the equalities 

Xi = psinai, yi = pcosai, (15) 



where < oti < tt/2. The five-dimensional volume element d£l p is 

dflp = sin 2 oti cos 2 Oj daj dxj dy^. (16) 

For the conjugate momenta, we proceed in a similar fashion introducing the hypermomentum k 
and the associated hyperangle /% 

ki = k sin Pi, qi = kcos/?j, (17) 

where owing to the energy conservation Q = k 2 j2\x. 

In the six-dimensional space, the kinetic energy operator T takes a separable form 

* 1 / d 2 5 d 1 r>2 \ 

where the hypermomentum operator ET is the generator of rotations in the six-dimensional space. 
The operator K 2 takes the form 

K 2 = -^ - 4cot(2a 4 )/- + -%- + -^-, (19) 

oaf oati sin aj cos z a« 

where the angular momentum operators i and Aj occurring in Q19|). are defined as 

k = -iXiX- — , Aj = -^X- — , (20) 

o Xi dyi 

and if and A 2 have eigenvalues £(£ + 1) and A(A + 1), respectively. The operator (fT9*|) has 
eigenvalues K{K + 4) where K = 2n + t + A for integer n. The quantum number K has been 
dubbed hypermomentum, its value is the same in all three Jacobi systems, and the corresponding 
eigenfunctions are known as hyper spherical harmonics (abbreviated HH hereafter). 

From now on we choose set 3 as the basic one and to simplify notation in the following we 
will suppress the label referring to our particular choice of the Jacobi system. The HH have the 
explicit form 

^KLM L 

where L is the total angular momentum resulting from vector addition L = I + A. In 1)21(1 
the symbol £l p denotes collectively all five angular variables with Yj m being the usual spherical 
harmonics. The square bracket in (|21[) indicates vector coupling of 1 and A producing a total 
angular momentum L and the appropriate quantum numbers associated with the latter operator 
are L and Ml. In 1)21(1 . the hyper angular eigenfunctions are 

x f (a) = Ng sin* a cos A aP^ +2 ' A+2 (cos 2a), 

where Pn (x) are the Jacobi polynomials (cf. [H]) and Nj^ is the normalization constant. The 
HH in ()21)) are orthonormalized using the volume element 1)16(1 and a stands for any of the a,. 
For clarity reasons, to avoid unnecessary complications, we are going to restrict our consid- 
erations to the simplest case when the particles have no internal degrees of freedom. When the 
particles are moving freely the spacial part of the three-particle wave-function will be described 
in the cm. frame by a plane wave, whose HH expansion can be written as 

e ,k.x + , q .y = (?*? J- y e AM L (n P ) 1 K V*PJk + 2{k P ) ^J^) 1 . (22) 



yKLM L («p) = *£(«) K(x) ® Y x (y)] L M T , ( 21 ) 



where Jk+2 denotes Bessel function of integer order and solid angle J7 K comprises five angular 
variables that specify the directions of the incident momenta in the six-dimensional space. It 
is worth noting that since the radial part in (122|) depends solely upon K, the plane wave is an 
invariant under six-dimensional rotations. The two-body interactions break up this invariance 
and for interacting particles the three-particle wave function in the continuum ^k jq (x, y) will 
have a more complicated HH expansion. For the simplest case of pair-wise central potentials 
depending upon the particle separations, we have 

* k , q (x,y) = |^E^W^p) *Wa(«,P) yl LuS^ ( 23 ) 

where the summation indices are L, Ml,K, £, A, K , I, A and the normalization condition reads 
^, q ,(x,y) t *k,q(x,y) d 3 xd 3 y = 5(k' - k) <J(q' - q). (24) 



/ 



In the expansion (|23|). only the total three-particle angular momentum L is a good quantum 
number and the as yet unspecified hyperradial part ipjL f * £-/t( k > p) must be determined from 
the underlying three-body dynamics. The dependence upon the tilded indices, associated with 
the incident momenta, enters solely via the asymptotic boundary condition at p — > oo. The full 
wave-function Vl/k,q(x, y) is a solution to the three-body Schrddinger equation 

(f + V r -Q)* k , q (x,y)=0, (25) 

with V = V\2 + V23 + V31 where V%j stands for the pair-wise interaction between particles i and 
j. Inserting Q23|) in (|25[l and projecting onto the hyperangular part of the wave function, results 
in an infinite set of coupled systems of differential equations enumerated by the conserved total 
angular momentum L - the only quantum number that does not mix. For an assigned L, we 
have 

p ' J (26) 



2p Y, {K'Z'\'\V\Kl\) i> L K , e , x , m (K,p), 



K'i'X' 



and we are left with a multi-channel situation where each channel is specified by three quantum 
numbers {Kl\}. Each system of equations (|26|) with L = 0, 1, 2, . . . is infinite because there is 
no upper limit for K, and, therefore, for practical reasons K must be truncated at some finite 
value K max so that the orbital momenta I and A are thereby restricted to vary in finite limits. 
The value of K max determines the order of the approximation and must be be large enough 
to ensure the convergence of the method. It should be noted that in (j26|) the potential term 
has been sandwiched between the HH functions y^LM (^p) anc ^ a fter integration over the five 
angular variables Q, p , these matrix elements depend solely upon the hyperadius p 

„ 3 3 

(K'£'\'\v\Ki\) = / y^' LMh {^ YH v ^p^p) y e KLM L ( n p) df V ( 27 ) 

J 3=1 i+i 

With the adopted here set 3 of HH, the computation of the potential matrix element V\i is 
relatively straightforward because the separation vector between particles 1 and 2 is proportional 
to |x| and the calculation of the potential matrix involves a single integration. By contrast, 
in the two remaining potentials the corresponding separations, r23 and r%i, respectively, are 



linear combinations of the two vectors x and y which unavoidably leads to five-dimensional 
integrations in the calculation of the potential matrix. Fortunately, there is a very efficient 
procedure to overcome this difficulty, based on the observation that all HH with fixed (K, L) 
values defined with respect to set i (cf. (|10[)) are linear combinations of HH belonging a set j 
and this transformation is effected by means of the so called Raynal-Revai (RR) coefficients [12| , 
viz. 

yKLM L ( a i'Xi>yi) = $^(*j \j\ii*i)KL ^JCLAf^, («i» *i> 9j), ( 2 8) 

where the RR coefficients (£j \j\li \)kl are functions of the angle specifying the kinematic 
rotation i — > j given in ()12|) . With the aid of (|28|) . the potential matrix of V23 may be computed 
in basis 1 where this task is simple and subsequently transformed to basis 3. 

The radial wave functions must be regular at p = and the boundary condition imposed on 
the solutions of (|26|) is 

^ L KexKix^P)-P K+5/2 for P^°- ( 29 ) 



For p — > 00 the potential term in IJ26JI goes to zero so that in this limit the system of equations 
((2l)|) becomes decoupled. In absence of the potential term, the asymptotic solutions of the 
radial equations are well known, they are given in terms of the Bessel and Neuman functions [Oj 
Jk+2(k>p) and Yk+2(k>p), respectively. With the Coulomb interaction present, for large p strong 
interaction potentials become negligible, and we set 

lim ^ L K£X m ( K , p) = 5 Kk 5 XX 8 d F c (r,, k P ) + (K£X\T L \K£X) H+(r,, K p), (30) 



with 



Ht(r),Kp) = Gc{r],np) + iF c (j],np). (31) 



where we have extended the standard definition of Coulomb wave functions so that the orbital 
quantum number £ which takes on integer values is replaced by C = K + 3/2 with half-integer 
C. The T-matrix occurring in (J3*U|) (K£X\T L \K£X) describes 3 — ► 3 scattering processes and 
is similar to that encountered in the two-body multichannel problem. In particular, it can be 
diagonalized which allows to determine the appropriate eigenphases 5lki\ whose rapid variation 
with energy indicates a resonance. 

3.2 Three-body final-state interaction enhancement factor 

The formalism developed in the preceding subsection will be now applied to calculate the en- 
hancement factor describing final-state interaction between three particles in a Borromean sys- 
tem. The approximation scheme will be based on the assumption that all interactions are of 
a short range and therefore the most important effect comes from the distortion of the wave 
function at very small separations. Bigger separation region is strongly suppressed by the small 
size of the interaction volume and gives relatively small contribution to the overlap integrals 
representing reaction transition amplitudes. Thus, similarly as in the two-body case, we shall 
define the enhancement factor as the limit reached by the square of the absolute value of the 
ratio - the full three-body wave function divided by the plane wave - when the size of the system 
represented by p goes to zero. It is apparent from (|26|l that the role of the centrifugal barrier 
plays now the quantum number K and the leading term in the wave function in the limit p — > 
has K = implying that £ = X = and L = 0. Therefore, we confine our attention solely to 
the L = set (|26|) where significant simplifications take place, namely we have X = £ and K 
necessarily even: K = 2(n + £). 



Formally, the three-particle enhancement factor will be defined as the limit 

2 



F(k, q) = lim 



*k, q (x,y) 



e «k-x+«q-y 

which by making use of (1221) and l)23j) . simplifies to the form 



(32) 



-F(fc, q) = lim 






2^ /*7il„(* n \ y Roo^ iK > 



K, 

oo en M2 



2 

(33) 



where we have dropped the uninteresting constant factor |3^joo(^p)l • ^-° obtain the enhancement 
factor it is sufficient to calculate the radial wave function ^ i?pp( K > P) f° r sman P an d take the 
limit indicated in (|33() . 

The calculation of the radial wave function is carried out similarly as for the two-body case 
except that scalar quantities now need to be replaced by matrices in channel space. These 
channels are specified by two quantum numbers K£ but to simplify the notation it is possible to 
use a single integer v = | (^K)(^K + 1) + £ + 1 enumerating all the states under consideration 
(^=1,2,3. . .). Truncating the infinite sequence of K values at some assigned K max , the total 
number of channels N is 

N=\{\K max + l){\K max + 2), (34) 

and we have to solve iV second order differential equations given in Q26JI . The physical solution 
of IJ2fi|) can be envisaged as a column vector in channel space. To obtain the physical solution 
we must first generate numerically N linearly independent particular solutions vanishing at the 
origin with an arbitrary slope at p = 0. These N solutions can be grouped together to form 
the N columns of a single N x N square matrix $ (hereafter we denote matrices by boldface 
symbols). In practice, this matrix may be obtained by solving N times the system of equations 
(|26|) imposing for small p the boundary condition 

®KtM^ P) = 6 kk S a M* +5/2 , (35) 

so that these particular solutions will be enumerated by two quantum numbers K£. Since the 
columns of $ are linearly independent, they span the space of all possible slopes. Therefore, 
any column vector solution ip of (|2fi|) must be expressible as some linear combination of the N 
columns of 3>, and we have 

il>{p) = *(p) • c, (36) 

with c denoting a column vector of constant coefficients. 

The matrix of physical solutions occurring in ()33|) differs from $ at most by a constant 
matrix and the latter will be explicitly determined by making use of the boundary condition for 
asymptotic p. For large p the physical solution must be of the form 

*(k,p) = F(k P )+H + (kp)-T(k), (37) 

where T is the true 3^3 scattering matrix. In formula l|37|) F denotes a diagonal matrix 
containing the regular solutions Fcij], up) of IJ26JI in absence of the strong interaction, whereas 
H + is another diagonal matrix containing outgoing hyperspherical waves Gc(i], Kp)+iF£,(j], np). 
Since $ is proportional to \I/, for some large p = p m , we may set the following matching condition 
for the wave functions and their derivatives 

*(/c,^»).C(k) = F(K Pm ) + H + (K Pm )-T(K), (38a) 

*(K,p m )'-C(«) = F( K p m y + H + ( K p m y-T( K ), (38b) 



where prime denotes derivative with respect to p and C(k) is a constant (albeit energy de- 
pendent) matrix to be determined. Eliminating T between the two matrix equations (|38|). we 
obtain 

C{k) = —[H+(K Pm ) ■ &(K, Pm y - H + ( KPm y ■ ^(K,p m )]-\ (39) 

and C(k) is seen to be proportional to the inverse of the Jost matrix. Clearly, formula (|39j) may 
be viewed as a three-body extension of (jHJ. Because the matrix C(k) involves the inverse of the 
Jost determinant it is bound to have the same singularity structure as the T-matrix. 
The physical solution can be now expressed entirely in terms of 3? 

*(k,p) = $(k,p)-C(k) (40) 

and the limiting procedure p — ► in (|33f) may be effected with the aid of (f3*5|) . Similarly as 
in the two-body case, the only terms in the summations giving non-vanishing contribution are 
those with K = 0. Leaving out irrelevant numerical factor, the enhancement factor in (|33|) takes 
the form 



F(k,q)=F(-k,-q) 



J2 c Ki( K )y e Koo( n ^ 



KC 



(41) 



where the expansion coefficients CKe(n) are provided as the first row of the matrix C(k) given 
in (HU). 

Formula (|41j) gives the enhancement factor in the form of an expansion in an orthonormal 
set of momentum dependent HH functions which turns out to be quite useful in effecting phase 
space integrations. We shall illustrate this point for the case when the transition matrix is 
assumed to be constant so that the cross section is essentially determined by the enhancement 
factor. The integration over the five angles df2 K , yields the total cross section 

a ^~TT-\ J>i«^)| 2 , (42) 

where n 4 results from phase space and f{n) denotes the incident flux factor. The summation 
over K extends over all even numbers from up to K max and in the following we take it to be 
the value at which the convergence has been attained. 

When the incident energy is fixed (i.e, k is a constant), the quantities of interest are usually 
the effective mass distributions, or equivalently, the corresponding kinetic energy distributions 
of the different pairs in their cm frame. Such distributions are obtained here by integrating 
only over the directions of k and q retaining the dependence upon the two remaining variables, 
k and (3. With our choice of 3 as the basic Jacobi set, the distribution of the center of mass 
energy T3 of the pair {12}, takes particularly simple form 

^r = VT3(Q-T 3 ) Y, WW CKtin) xUP) x£(/3), (43) 

3 K'Kl 

where T3 = Q sin and the square root factor is a remanet of the phase space. All other 
factors that depend only upon n have been dropped as they may be absorbed in the arbitrary 
normalization constant. The calculation of cm kinetic energy distributions of the two remaining 
pairs, i.e {23} and {31}, respectively, may be carried out in exactly the same manner but the HH 
occurring in (|4*Tj) need to be first transformed to the appropriate Jacobi system by applying a 
kinematical rotation. Thus, using the transformation (|28|1 in (j41j) . the distribution of the kinetic 
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energy of the {23} pair, T\, is 

ai1 K'Ki' 3 i 3 h ' (44) 

x {£i£i\£s£3)ko (t 3 £z\h£i)K'o, 

where T\ = Q sin (3\ and the corresponding distribution of T2 follows from 1)44(1 in result of 
Pi 2 permutation. It should be perhaps clarified that this permutation in general changes the 
kinematic rotation angle that enters the RR coefficients and in consequence the T2 distribution 
does not have to be the same as that given in (|44l) . 

4 Comparison with experiment and conclusions 

In the preceding section our considerations have been quite general and now we wish to apply 
this theory to the investigation of FSI effects in the pp — ► ppr] reaction close to threshold. 
To simplify matters we assume that the excitation energy is sufficiently low so that the total 
orbital momentum in the final state is zero and the 3 Po — > 1 Sqs transition amplitude is the 
sole contributor to the cross section. We label the two protons as 1 and 2, respectively, and 
the meson as 3 choosing Jacobi set 3 for all computations. With the two protons in a singlet 
state, Pauli principle requires that in this frame the orbital momenta £3, A3 take only even 
values. This guarantees that the spatial part wave function will be symmetric with respect 
to the {12} permutation. We have tried several possible forms of the pp potential in the 1 Sq 
state: a delta-shell, a Gaussian or a fully realistic Reid potential (for details cf. [5]). The 77-p 
interaction is poorly known, for the real part of the scattering length values between 0.2 and 
about 1.0 fm have been suggested ^3]. Additional difficulty stems from the multi-channel nature 
of this interaction: already at threshold there are open channels so that r/p scattering length 
is a complex number. At present there is not enough information to include these additional 
channels into our formalism, therefore in this work rj-p interaction has been simulated by the 
simplest non-absorptive delta-shell potential operative in the s-wave only whose range has been 
arbitrarily fixed to be 1 fm. The depth of this potential can be then adjusted as to yield an 
assigned value of the r/p scattering length a vp . Since the precise value of the latter is not available 
|13j we used three values 0.5, 1.0 and 1.5 fm, respectively, which we believe is a representative 
sample. 

All our computations have been carried out for the lowest excitation energy equal Q=15.5 
MeV for which two-particle invariant mass spectra are experimentally available. The system of 
equations (|2T)|) was perpetually solved by increasing successively in each step the value of K max 
until convergence has been attained. For delta-shell and Gaussian pp potentials this occurred 
at K max = 24 but with Reid potential this figure must be doubled. Similarly as in the two-body 
case jH] the results are completely insensitive to the shape of the pp potential. The results of 
our computations are presented in Fig. |2] where they are compared with the data from |3j. Not 
unexpectedly, the dominance of the pp interaction is apparent from both plots. Since the values 
of a vp > 1 fm are probably rather unrealistic |13j , we may say that the r]p interaction is of 
marginal importance as far as invariant mass distributions are concerned. 

Comparing the invariant mass spectra obtained in the two-body approach (Fig. ^) with 
those resulting from the full three-body calculation (Fig. |2J) we can see that even when the r/p 
interaction is completely disregarded, the invariant mass plots are different in these two cases. 
Although the input in both approaches is the same but the underlying calculational schemes are 
different. Owing to the proper boundary condition l|37() also in absence of ij-p interaction the 
three-body wave function has entangled form, i.e it cannot be expressed as a product of the p-p 
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Figure 2: Invariant mass distributions at Q=15.5 MeV obtained from a three-body calculation: left 
panel pp pair; right panel r/p pair. The data are from 0|. 



wave function times a plane wave associated with the free propagation of the r\. Since in the 
approximation using the two-body enhancement factor the very existence of the r\ is unaccounted 
for, the latter does not depend upon the meson kinematics (in [5] we made an attempt to lift this 
deficiency by introducing ad hoc a linear dependence upon q). By contrast, even in absence of 
r?-p interaction the three-body enhancement factor accounts for the 7/ propagation. In particular, 
the mass of the r\ strongly influences the dynamics of the three-body system as can be seen from 
(HSU and (1271) . 

The three-body calculation is in both spectra closer to experiment. We wish to note here 
that the curves presented in Fig. [2] and as well as the dashed curves from Fig. ^ contain no 
adjustable parameters except for the overall normalization. Our three-body calculation clearly 
favors smaller a^p values a np ~ 0-r0.5 fm. For a vp > 0.5 fm the height of the close to threshold 
peak in s pp plot in Fig. |2] is depressed which is accompanied by a build up of a pronounced 
shoulder at the high-energy end. At the same time another shoulder appears at the low-energy 



end in s^p distribution in Fig. |2 All these features worsen the agreement with experiment. 
Unfortunately, it would be difficult to improve the existing estimates of the o w scattering length 
using the data displayed in Fig. especially that absorptive effects have been ignored. 

Summarizing, we have developed a general three-body formalism for calculating FSI effects 
in a three-particle final state. To the best of our knowledge, this is the first attempt in the 
literature to derive a three-body final-state interaction enhancement factor. The presented 
calculational framework employs the hypersherical harmonics method and is applicable for three- 
body systems in which no binary bound state can exist. The necessary input which must be 
provided are all three pairwise potentials. Detailed computations carried out for the pp — ► ppr\ 
reaction at Q=15.5 MeV show that in the invariant mass spectra the role of the 77-p interaction 
is marginal and more important is the proper boundary condition imposed on the final-state 
wave function. The three-body calculations involving no adjustable parameters reproduce quite 
well the experimental invariant mass spectra. 

Partial support under grant KBN 5B 03B04521 is gratefully acknowledged. 
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